Exact solution of the nuclear pairing problem 
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In many applications to finite Fermi-systems, the pairing problem has to 
be treated exactly. We suggest a numerical method of exact solution based 
Q I on SU(2) quasispin algebras and demonstrate its simplicity and practicality. 

■ We show that the treatment of binding energies with the use of the exact 

pairing and uncorrelated monopole contribution of other residual interactions 



(N 



\ can serve as an effective alternative to the full shell-model diagonalization in 

pLn ' spherical nuclei. A self-consistent combination of the exactly treated pairing 

Q ■ and Hartree-Fock method is discussed. Results for Sn isotopes indicate a good 

CN , agreement with experimental data. 

> 

^ \ Pairing correlations play an essential role in nuclear structure properties including bind- 

^ ■ 

ing energies, odd-even effects, single-particle occupancies, excitation spectra, electromag- 

O 

^ \ netic and beta-decay probabilities, transfer reaction amplitudes, low-lying collective modes, 
level densities, and moments of inertia |[l|-p[. The revival of interest in pairing correlations 



^ \ is related to studies of nuclei far from stability and predictions of exotic pairing modes 



Metallic clusters, organic molecules and FuUerenes are other examples of finite Fermi 



^ ' systems with possibilities for pairing correlations of the superconducting type 



The conventional description of pairing usually employs the classical BCS approach 
used in theory of superconductivity. This approximate solution has a very good accuracy for 
large systems and becomes exact in the asymptotic limit The shortcomings of the BCS 
approximation for small systems are well known, see for example and references therein. 
The major drawback of the BCS is the violation of particle number conservation, which gives 
rise to deviations from the exact solution for small systems. Various ideas have been sug- 
gested to correct this deficiency, such as the number projection mean-field methods p0|-p!2 



coherent state approach [|T^, stochastic number projection |T^, statistical descriptions []T3 



treatments of residual parts of the Hamiltonian in the random phase approximation ||T^,|T7 



1 



and recurrence relation methods [T^JT^. These methods have found only a limited num- 
ber of practical applications; for some approaches the obtained results did not manifest the 
desired accuracy whereas other methods are limited by practical complications. BCS-like 
approximate theories have a number of other deficiencies when applied to small systems . 
In particular, in the region of weak pairing the BCS has a sharp phase transition from the 
paired condensate to the normal state with no pairing (trivial or zero gap solution), whereas 
exact solutions exhibit the existence of exponentially decreasing pairing correlations all the 
way down to the zero pairing strength. This difficulty makes the BCS method unreliable 
for applications to weakly bound nuclei and calls for improvements and extensions such as 
BCS+RPA py]. There are also serious problems related to the correct description of pairing 



in excited states. 

The exact pairing (EP) method presented in this work allows one to solve exactly the 
general pairing Hamiltonian 

H ^ y ^i^i m^j m ~l~ ^ ^ j j' ^ y (^jm^jm^j'm'^j'm' ; ^jm — ( 1) —m j (1) 

jm j^j' m,m' 

where {e^} is the set of single-particle energies, diagonal Gjj are pairing energies, and 
Gjj/ = Gj> j for j 7^ j' are pair transfer matrix elements, {j' j). The practical usefulness 
of the EP algorithm comes from the facts that it is exact, fast and allows a straightforward 
extension for an approximate treatment of other components of residual interactions. Some 
realistic examples are presented in this work to emphasize these points. As a result, it 
becomes unnecessary to use complex approximate methods, such as BCS and all those 
associated with it, when exact results that are free of all problems discussed above can be 
obtained with an almost equal or even smaller effort. 

A number of methods for treating the pairing problem exactly have been previously 
proposed. The Richardson method, described in the series of papers ||21|, provides a for- 
mally exact way for solving the pairing Hamiltonian. This method reduces the large-scale 
diagonalization of a many-body Hamiltonian in a truncated Hilbert space to a set of cou- 
pled equations with a dimension equal to the number of valence particles. Recently, exact 



solutions have been approached by introducing sophisticated mathematical tools such as 
infinite-dimensional algebras ||22[ . Such formally exact solutions have a certain merit from a 
mathematical point of view and for developing and understanding approximate calculations. 
However, due to their complexity they are not very useful in solving practical problems in 
nuclear physics. 

The natural way of solving the pairing problem is related to the direct Fock-space di- 
agonalization. For deformed nuclei with the doubly degenerate single-particle orbitals this 
approach supplemented by the appropriate use of symmetries and truncations was already 
shown to be quite effective |2^j2^. Our goal is to combine the exact treatment of pairing 



with the approximate inclusion of other parts of the residual interaction. Here the framework 
of the rotationally invariant shell model is the most convenient, especially because it allows 
us to fully utilize the well known ideas (see for example based on the existence 



of quasispin symmetry in paired systems first studied in the 1940's by Racah [Id]. In the 



context of the shell model similar ideas were utilized in |]2^. We use the quasispin algebra 
for each subset of degenerate single-particle levels. The fact that Racah's degenerate model 
is analytically solvable comes purely from this algebraic feature. The spherical shell model 
with its m-degeneracies is a perfect arena for applying this method and therefore we will use 
notations associated with the case of spherical symmetry and j-j coupling. This certainly 
does not limit the generality of approach. 

The Hamiltonian Eq. (|ID can be rewritten as 

H = Y: e,n, + 2 Y: + E G^rLjlr, , (2) 

3 3 3 j' 

by introducing the partial quasispin operators Lj', Lj and Lj for each j-level as follows 

^ ra ^ ra 

^ = \Y. i^lm^^rn = ^N, - , (4) 



2 ^ V ' 2J 2 
where Nj is the particle number operator and Qj = (2j + l)/2 is the pair degeneracy of a 



given single-particle level j . It can be shown directly from the definitions of L j , Lj and Lj 
that they form an SU(2) algebra of angular momentum, 

Expressed in terms of quasispins, the Hamiltonian in makes the pairing problem equiva- 
lent to the problem of interacting spins in a magnetic field, a generalized form of the Zeeman 
effect [^. It is clear that every square of the partial quasispin = Lj'Lj — Lj + {L 



z\2 



commutes with the Hamiltonian making Lj, corresponding to the eigenvalue Lj{Lj + 1) of 
L^, a good quantum number in the pairing problem. This brings in the major simplification 
of the problem. The maximum value that Lj can take is Qj/2 , which happens for a fully 
paired subshell, such as, for example, the completely occupied subshell case where A'^ = 2Qj, 
Lj = Lj = Lower values of the quasispin quantum number correspond to the Pauli 

blocking of a part of the pair space Qj by Sj unpaired particles. This reduces the allowed 
space to Qj — sj . The number sj can be called the seniority of a given j-shell. Being related 
to the quasispin, sj are conserved by the pairing interaction (0). Introducing the partial 
seniority Sj and partial occupancy Nj by means of 

L, = l{n,-s,) L^^ = ^{N,-n), (6) 

we will use \sj,Nj) instead of the SU(2) notation \Lj,Lj) . 

The quasispin projections Lj , or the partial occupancies Nj , are not conserved because 
of the pair transfer term Gjj> (pair vibration). The usual constraints on angular momentum 
^ \Lj\ < Lj lead to Sj < Nj < 2Qj — Sj and Sj < Vtj that have obvious interpretations. 
Furthermore, since both Lj and L| must be simultaneously either integers or half-integers, 
Sj is of the same parity as Nj , confirming that particles are transfered only in pairs. Finally, 
there is a given total number of particles in the system = Y^j Nj , and one can introduce 
the total seniority s = J2j Sj . Both quantities are conserved by the Hamiltonian and must 
be of the same parity. 

For each representation given by the set of quantum numbers {sj} , one can construct 
basis states {Nj}) going through all permutations of fermions allowed by above 



constraints. Finally, a Hamiltonian matrix can be constructed in this basis using standard 
properties of the angular momentum operators 

L±|L, L,) = ^{LtL,){L±L, + 1) |L, L ± 1) . (7) 

Diagonal elements become 

({5,}, {N,}\H\{s,}, {iV,}) = + ^{N, s,)(2n, - s, - + 2)) , (8) 

3 

and off-diagonal elements that transfer pairs are 

{{s,},...N,+2,...N^,-2,...\H\{s,},...N,,---N,,,...) 

= - ^^■')(2^" - V - N,, + 2)(2n, - s, Nj){N, s, + 2) . (9) 

Diagonalization of this matrix for each representation, a given set of partial seniorities {sj} , 
is the final step in the solution. The largest Hamiltonian matrix in partial seniority basis, 
which corresponds to the lowest allowed total seniority, has a dimension much lower than 
the total fermionic many-body space. Even for the set of valence orbits encountered in 
heavy nuclei it does not exceed several thousands. As seen from the above discussion, the 
Hamiltonian matrix is very sparse. 

Each state of a non-zero seniority is degenerate since Sj unpaired particles are un- 
touched by the Hamiltonian and are free to move within a given j-shell, provided that 
they all remain unpaired. The number of these fermionic degrees of freedom for each j is 
C(2Qj, Sj) — C(2Qj, Sj — 2) where C(m, n) — m\/{n\{m — n)\) is a binomial coefficient. The 
total degeneracy of each state in the {sj} seniority class is 

l[[C{2n„Sj)-C{2n„s,-2)] . (10) 

3 

The resulting degenerate states that carry seniority quantum numbers can be further classi- 
fied according to other symmetry groups of the Hamiltonian such as angular momentum. All 
total seniority zero states have spin zero. For non-zero seniorities, one first has to determine 
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allowed angular momenta for a given subshell j with sj unpaired particles and then make 
all possible couplings of subshells to a total spin. 

To illustrate the practical application of this algorithm, we show below examples involv- 
ing chains of Ca and Sn isotopes. The Ca isotopes occupy the fp shell with the model space 
consisting of four levels fj/2, P3/2, f5/2, and pi/2 ■ This amounts to a total neutron capacity of 



20. The pair transfer matrix elements are taken from the FPD6 interaction |2^. The results 
in Fig. |l|(b) are obtained with single-particle energies appropriate for the above-mentioned 
j-levels in ^^Ca: -9.9, -5.1, -1.6, and -3.1 MeV, respectively. Fig. |l|(b) shows the correlation 
energies in the Ca isotopes, defined as 

i5„„=£;-i:.,iv,-i:-^M|^, (11) 

where E is the energy of the ground state and Nj is the (noninteger) expectation value 
of the ground state occupancy of the j-level, found as a result of the EP calculation. It 
is known [|l| that in a system with a discrete single-particle spectrum the nontrivial BCS 
solution exists only at a sufficiently strong pairing interaction which can overcome the single- 
particle level spacings. In contrast to that, in macroscopic Fermi-systems the nontrivial 
Cooper phenomenon exists at any strength of the attractive pairing interaction. The fact of 
inadequacy of the BCS approximation for weak pairing was pointed out by many authors, for 
example P,p^. Near the ^^Ca shell closure there is a 4.8 MeV gap between /7/2 and higher 
single-particle orbitals, and as it can be seen from Fig. |l](b), BCS can no longer support the 
pairing condensate. It gives instead a normal Fermi-gas solution with zero pairing energy, 
whereas in reality the EP solution demonstrates significant pairing effects with almost 2 
MeV condensation energy. 
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FIG. 1. Panel (a): neutron separation energies in Ca isotopes, EP+monopole treatment of other 
interactions is compared with experimental points (diamonds); panel (b): pairing correlation energy 
in even-even Ca isotopes, the exact EP calculation (solid line) is compared with the standard BCS 
(dashed line). 

In most of the real cases the pure pairing interaction ( J = 0, T = 1 ) is strongly modified 
by residual interactions in other channels. A significant part of the remaining correlations 
in spherical nuclei can be treated by including only monopole (diagonal) part 

j 3+3' 

where 

y33' = ,0/90^"' A + 1)^^' ^'^"i^i^' ^■^■') (13) 

is the monopole-monopole part of the total interaction V written in eq. (|13]) in terms of 
the diagonal matrix elements for the pairs Pairing (J = 0; j j\V\J = 0; j' j') = 

G j ji ^JVtjVtji is not included in the sum since it is treated exactly by EP and its monopole 
term is present in Eq. (^). One- neutron separation energies shown in Fig. ^(a) were calcu- 
lated with this EP plus monopole method, using the FPD6 interaction and single-particle 
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energies for '^^Ca. These results agree very well with experimentally observed separation 
energies, shown on the figure with diamonds. Ground state energies, obtained with this 
method for Ca isotopes differ from those from the exact shell model diagonalization, using 



the total FPD6 interaction, by less than 0.5 MeV; certainly the EP plus monopole 
results are exact for full and empty shells and for one-particle or one-hole cases. 

For the second example of application we show in Fig. ^(b) the results for the Sn isotopes. 
The model space here consists of five single-particle levels /iii/2, (^3/2, -51/2, g7/2, and ^5/2, 
which can accommodate up to 32 neutrons. We will be interested in the evolution of neutron 
separation and correlation energies along the chain of the isotpes rather than in the specific 
spectra which were studied earlier with a number of approximate methods, for example ||31|| , 



and can be considered in the future with our exact approach. The pairing matrix elements 



were obtained from the G-matrix derived from the recent CD-Bonn |^ nucleon-nucleon 
interaction with ^^^Sn as a closed shell, where the Q-box method includes all non- folded 
diagrams up to the third order in the interaction and sums up the folded diagrams to 
infinite order The most complex case is ^^^Sn (half-filled shell) with 601,080,390 many- 



body states of which 272,828 are of spin zero, so that even with the use of the angular 
momentum projection the problem remains difficult for direct diagonalization. There are 
only 420 independent seniority sets in the seniority basis. The largest matrix, s = , has a 
dimension of 110 , the diagonalization of which is a trivial problem. Solution of the pairing 
problem for ^^^Sn with the algorithm discussed above is therefore simple and extremely fast. 

The correlation energies obtained with the above G-matrix and single- hole energies -9.76, 
-8.98, -7.33, -7.66, and -7.57 MeV, based on the data for the ^^^Sn isotope pO[, are shown 



in Fig. 0(b). A reduction of pairing also happens between gj/2, and the rest of the 
single-particle orbitals /iii/2, (^3/25 Si/2 , but here BCS is just weakened which results only in 
relatively small suppression of correlation energy. A related application of EP is shown in 
Fig. 0(a) where one-neutron separation energies are calculated for the Sn isotopes, including 
those beyond ^^^Sn. These energies were obtained from the fully self-consistent spherically 
symmetric solution of Hartree-Fock (HF) equations, using the SKX interaction ||3^, with 
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the EP solution based on the above G-matrix at each HF iteration. The HF+EP method 
works as follows. An initial guess for the HF potential is used to obtain a set of spherical 
single-particle energies and densities for all of the occupied and valence orbitals. The single- 
particle enegies for the valence space plus a fixed set of two-body matrix elements are used 
in the EP method to obtain the pairing correlation energy and the valence single-particle 
occupation numbers. The initial set of single-particle densities plus the pairing occupation 
factors are used to obtain a new potential with the Skyrme hamiltonian. This proceedure 
is iterated until convergence (typically 60 iterations). The total energy is the sum of the 
Skyrme HF energy and the pairing correlation energy Eq. (0). The plot H(a) exhibits 
an odd-even staggering that can only be attributed to pairing. The agreement between 
experiment and theory is excellent given that no parameters have been adjusted for these 
particular data. 
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FIG. 2. Panel (a): neutron separation energies calculated for the Sn isotopes using the 
self-consistent solution of Hartree-Fock plus Exact Pairing. Experimental points are shown for 
comparison. Panel (b): pairing correlation energy in even-even Sn isotopes, EP (solid line) and 
BCS (dashed line). 

The EP algorithm presented here in our view has a good future as a tool for different 
calculations related to pairing. It is exact, fast and reliable which makes it perfect for pure 
shell model calculations with pairing, fast estimates of binding energies and spectroscopic 
factors, for the use as a basis for treatment of other interactions, iterative HF calculations 
and many other tasks. 
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